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Abstract 

Utilizing recent DIS measurements (i*2,L) and data on dilepton and high-£/r jet 
production we determine the dynamical parton distributions of the nucleon gener- 
ated radiatively from valence-like positive input distributions at optimally chosen 
low resolution scales. These are compared with 'standard' distributions generated 
from positive input distributions at some fixed and higher resolution scale. It is 
shown that up to the next to leading order NLO(MS, DIS) of perturbative QCD 
considered in this paper, the uncertainties of the dynamical distributions are, as 
expected, smaller than those of their standard counterparts. This holds true in 
particular in the presently unexplored extremely small-x region relevant for evalu- 
ating ultrahigh energy cross sections in astrophysical applications. It is noted that 
our new dynamical distributions are compatible, within the presently determined 
uncertainties, with previously determined dynamical parton distributions. 



1 Introduction 



The parton distributions of the nucleoli are extracted from deep inelastic scattering data 
by essentially two different approaches which differ in their choice of the input distributions 
at some low scale Qq. In the common approach, e.g. [TJ [2j [3] , Qq is fixed at some arbitrarily 
chosen Qq > 1 GeV and the corresponding input distributions are unrestricted, allowing 
even for negative gluon distributions [31 HJ E] in the small Bjorken-x region, i.e. negative 
cross sections like F^(x, Q 2 ). Alternatively [6j [7J |8] the parton distributions at Q £ 1 
GeV are QCD radiatively generated from valence-like positive input distributions at an 
optimally determined Qq = \i < 1 GeV (where 'valence-like' refers to af > for all input 
distributions xf(x, /x 2 ) oc x af (l— x) bf ). This more restrictive ansatz implies, of course, less 
uncertainties concerning the behavior of the parton distributions in the small-x region 
at Q > fi which is entirely due to QCD dynamics. In particular it provided unique 
(steep) predictions [HI IS] for the experimentally unexplored region x < 10~ 2 which were 
subsequently first confirmed in [T0| [TTj . 

This predictive power is especially important for investigations concerning cross sec- 
tions [12J of ultrahigh energy particles (neutrinos) produced via astrophysical acceleration 
processes, e.g. in active galactic nuclei, black holes or in the decays of very massive par- 
ticles (see, for example, [THJ [HI US])- Here one needs a somewhat reliable knowledge of 
parton distributions at the weak scale Q 2 = down to x ~ 1CT 9 (x ~ jlm^Ey) at 
highest energies of E v ~ 10 12 GeV which requires extrapolations into the yet unmeasured 
small-x region x < 10~ 3 . Furthermore this 'radiative' approach based on QCD dynam- 
ics is also useful for connecting nonperturbative models valid at Q < 1 GeV (like chiral 
quark-soliton models [16-20] and statistical parton models [21-24]) with the actually 
measured distributions at Q > 1 GeV. 

In the present paper we confront recent precision deep inelastic scattering (DIS) ep 
data, Drell-Yan dilepton and high-p T jet data with radiatively generated parton distri- 
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butions arising from a valence-like positive input at Q < 1 GeV, following and extending 
the latest GRV98 analysis [8]. Moreover we study the dependence and stability of the 
small-a; predictions, in particular of their extrapolations down to 1CT 9 , with respect to 
a different choice of the factorization scheme (MS versus DIS). Furthermore we compare 
these 'dynamical' results with the ones obtained from the common evolution approach 
being based on a non-valence-like input at Qq > 1 GeV. In particular we shall compare 
their associated uncertainties. As should be clear by now, it will turn out that these uncer- 
tainties are indeed smaller for the radiatively generated parton distributions, particularly 
in the small- x region, due to their valence-like input and the sizeably larger evolution 
distance starting at Q < 1 GeV. 



2 Formalism 



The aforementioned analyses are undertaken at the next-to-leading order (NLO) of per- 
turbative QCD within the modified minimal substraction (MS) factorization and renor- 
malization scheme. For the radiative model we shall also present results as obtained 
within the so-called DIS factorization scheme [221 EE] • Heavy quarks (c, b, t) will not 
be considered as partons, i.e. the number of active flavors rif appearing in the splitting 
functions and the corresponding Wilson coefficients will be fixed, n/ = 3. This defines 
the so-called 'fixed-flavor number scheme' (FFNS). As argued in [27], it is nevertheless 
consistent and correct to utilize the standard variable rif scheme for the (3— function, and 
we shall adopt this procedure in our evaluation of the running coupling constant a s {Q 2 ). 
Up to NLO, the strong coupling a(Q 2 ) = a s {Q 2 ) / An evolves according to 

da/dlnQ 2 = -f3 a 2 - /^a 3 (1) 

where /3q — 11 — 2n//3 and (5\ = 102 — 38n//3. Here we utilize the exact numerical 
(iterative) solution for a(Q 2 ) since it is mandatory in the low Q 2 region [H] relevant for 
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the valence-like approach. The exact solution of (1) can be written implicitly 

Since /3 ,i are not continuous for different n/, the continuity of a{Q 2 ) requires to choose 
different values for the integration constant A for different flavor numbers rif, A^ n f\ which 
are fixed by the a(Q 2 ) matchings at Q = m Cy b,t- We have chosen 

m c = 1.3GeV, m b = 4.2 GeV, m t = 175GeV, (3) 

which turn out to be the optimal choices for all our subsequent LO and NLO analyses 
of heavy quark production. This exact solution reduces to the common approximate 
'asymptotic' solution 

W J " A,ln(Q 2 /A 2 ) (31 [ln(QVA2)]2 ^ > 

which turns out to be sufficiently accurate for Q 2 >, 2 GeV 2 and which, moreover, is easier 
to use for practical applications of our results. The values for Aw', as well as for A^ n f\ 
corresponding to our various LO and NLO global dynamical fits are given as follows: 
In LO, where ft = and = A^\ A { ^ 6) = 210.3, 181.8, 138.3, 70.1 MeV; in 

NLO(MS) AS^' 5 ' 6) = 302.8, 251.0, 172.8, 70.0 MeV and A^ 5,6) = 269.7, 184.5, 72.9 MeV; 
in NLO(DIS) ASs' 5 ' 6) = 288.5, 238.0, 162.9, 65.6 MeV while A^g' 6) = 255.9, 173.9, 68.3 
MeV. 

Let us now turn to the update of our LO and NLO (MS) GRV98 distributions [8] 
which consists of a fine tuning of the valence-like input densities xf(x, Qq) as well as of 
the input scale Q = \i < 1 GeV. The non-singlet input densities u v , d v , A = d — u and 
the valence-like input distributions d + u, s = s and g in the singlet sector are generically 
parametrized as 

xf(x,Ql) = N f x a f(l-x) b f(l + A f ^c + B f x) (5) 
subject to the constraints u v dx = 2, J d v dx = 1 and 

/ x[u v + d v + 2(u + d + s) + g] dx = 1 . (6) 
Jo 
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Since the data sets we are using are insensitive to the specific choice of the strange 
quark distributions, we continue to generate the strange densities entirely radiatively [S] 
starting from s(x,Qq) = s(x,Qq) = in the valence-like approach where Qq < 1 GeV. 
In the common standard approach where Qo > 1 GeV we choose as usual s(x,Qq) = 
s(x,Qq) = [u(x,Qq) + d(x,Qo)]/A. Furthermore, since all our fits did not require the 
additional polynomial in (5) for the gluon distribution, we have set A g = B g = 0. This 
left us with a total of 21 independent fit parameters, including a s . 

These free parameters have been fixed using the following data sets. The HERA ep 
measurements [28j EH] for Q 2 > 2 GeV 2 for the 'reduced' DIS one-photon exchange cross 
section a r = F 2 — (y 2 /Y + )F L together with the full neutral current (7, Z°) sector [3U] 

(j e±p (r O 2 ) = ( 2na2Y + \ d2(T NC _ f nc _ f_ f nc Xz. T pNC (7) 
where Y± — 1 ± (1 — y) 2 and 

^2 N l = F 2>L -v e KF2l + (v 2 e + al)K 2 F 2 z h 

F 3 NC = -a e nFf + (2v e a e )n 2 Fi (8) 

with v e = —\ + 2sm 2 6 w , a e = = 4 sin 2 9 W cos 2 9 W (Q 2 + M%)/Q 2 using 

sin 2 ^ = 0.2312 and Mz = 91.1876 GeV. Note that the structure functions in (7) re- 
fer to the radiatively corrected ones as presented by the experimentalists. Furthermore, 
the well-known standard target mass corrections to F 2 have been taken into account in 
the medium to large x-region for Q 2 < 100 GeV 2 . Since the experimental extraction of 
the usual (one-photon exchange) F2(x,Q 2 ) from d 2 a/dx dQ 2 is obviously (parton) model 
dependent, we have chosen to work with the full NC framework in order to avoid any 
further dependence on model assumption. However, it turned out that fitting just to 
F 2 (x,Q 2 ) gives very similar results. Furthermore, since the longitudinal structure func- 
tion Fl contributes sizeably only for large inelasticity y, in most of the kinematic range 
we have a r ~ F 2 . In addition we have used the fixed target F 2 data of SLAC [31], 
BCDMS [32J , E665 [33] and NMC [34] subject to the standard cuts Q 2 > 4 GeV 2 and 
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W 2 = Q 2 {\ — 1) + rrip > 10 GeV 2 , together with the structure function ratios F^jFl 
of BCDMS [35J, E665 [3S] and NMC [37]. The data for heavy quark (c, h) production, 
being theoretically described in the fixed-flavor number factorization scheme by the fully 
predictive fixed-order (LO/NLO) perturbation theory to be discussed below, are taken 
from [3"8l EHJ HO] for F% and from [3D] for F%. Furthermore the Drell-Yan muon pair pro- 
duction data of E866/NuSea [H] for d 2 a pN /dx F dM^+^- with N = p,d have been used, 
as well as their asymmetry measurements [42 J for a pd /a pp . These data are instrumental 
in fixing d — u (or d/u). The relevant LO/NLO differential Drell-Yan cross sections can 
be found in the Appendix of [IB], except for eq. (A8) which has to be modified [HJ H5] 
in order to conform with the usual MS convention for the number of gluon polarization 
states 2(1 — e) in 4 — 2e dimensions. Finally, the Tevatron high-px (or Et) inclusive jet 
data of DO [IB] and CDF [IT] have been used together with the fastNLO [IE] package for 
calculating the relevant cross sections at NLO. It should be mentioned that all these data 
sets correspond to a total of of 1739 data points. 

As already noted, the LO and NLO heavy quark contributions Ff ,b are calculated in the 

FFNS and contribute to the total structure functions as Fi(x, Q 2 ) = E^ ght + i^ heavy where 

'light' refers to the common u, d, s (anti)quarks and gluon initiated contributions [44J, and 

pheavy _ pc _|_ pb Top quark contributions are negligible. The LO 0(cn a ) contributions^] 

to F 2 h L , due to the subprocess 7*g — > hh, have been summarized in [7J, and the NLO 

0(a 2 ) ones are given in [5T]|52]. These contributions are gluon g(x,n 2 F ) dominated where 

the factorization scale should preferably be chosen [53] to be [i 2 F = 4m 2 . As we shall 

see the resulting predictions are in perfect agreement with all available DIS data on 

heavy quark production and are futhermore perturbatively stable [53]. Even choosing 

a very large scale like fj, F = 4(Q 2 + 4m 2 ) leaves the NLO results essentially unchanged 

[8| 154"] in particular at small-x. This stability renders attempts to resum supposedly 'large 

1 It has become common to consider the 0(a s ) contributions to F^ sht as LO [35] . Therefore the 2-loop 
0(a 2 s ) Wilson coefficients [33] are required for i^ lght at NLO. (Simplified parametrizations of the relevant 
Wilson coefficients can be found in [50].) 
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logarithms' (In Q 2 /ml) in heavy quark production cross sections superfluous. The heavy 
quark contributions to the remaining NC structure functions and _F 2 Z L in (5) are 
quantitatively negligible. Nevertheless, they can be simply obtained from the 7-exchange 
contributions F 2L by substitutions like e 2 — > 2e q v q and e 2 — > v 2 + a 2 , respectively, where 
v q = ±| — 2e g sin 2 W and a q = ±| with ± referring to an up- or down-type quark. The 
heavy quark contributions to F^ z,z vanish in LO [55| 156] , and are negligibly small in NLO 
since effectively F% ~ h — h at the relevant large values of Q 2 . 

In order to test the dependence of our results on the specific choice of the factorization 
scheme in NLO, other than the commonly used MS scheme, we also perform our NLO 
analysis using the deep inelastic scattering (DIS) factorization scheme [231 |2"6"] . Here the 
MS Wilson coeflients are absorbed into the parton distributions, or more precisely into 
their evolutions, i.e., into the splitting functions. This transformation to the DIS scheme 
in NLO is achieved via [H] 

p(l) __. pW _ p(l) , a Ar W (q\ 

r NS "'NS.DIS ~~ -'NS 1 A J 0^^2,NS V*) 



(10) 



where 



ACSs = -<$. • *5« = - f % • % I ■ (ID 

The light u, d, s quark contributions to F$, for example, in the NLO (DIS) factorization 
scheme now simply become 

F^ ht (x,Q 2 )=x £ e 2 q [q(x,Q 2 ) + q(x,Q 2 )] ms . (12) 

q=u,d,s 



The quantitative difference between the NLO (MS) and NLO (DIS) results will turn out 
to be rather small. Having obtained the parton distributions ^q\x, Q 2 )bis an d g{x, Q 2 )ms 
from an explicit NLO analysis of F 2 (x,Q 2 ) in the DIS factorization scheme, one can 
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transform them to the MS scheme via (see [7J, for example) 



{ q\x,Q 2 ) 



g(x,Q 2 



V(:r, Q 2 )dis - a 



(x,Q 2 ) + 0(a 2 ) (13) 



s(x,Q 2 ) D is + a c2®S DI s + ^®g DIS (x,Q 2 ) + 0(a 2 ) 



(i) 



where 



c£>(*) 



i + -: 2 / |n 1 - : 



4n ; 



^ + (1 



In 



+ -(9 + 5z) 



- 1 + 8z(l 



(14) 



(15) 
(16) 



with rif = 3. This transformation to the MS scheme then allows also for a consistent NLO 
analysis of heavy quark and Drell-Yan dimuon production processes in the DIS scheme, 
using their well known theoretical MS expressions, as well as for a consistent comparison 
of our DIS results with the ones obtained in the MS factorization scheme. 



2a. Estimates of uncertainties 

Our evaluation of the parton distribution uncertainties is based on the Hessian method 
with the Hessian matrix defined via 

AX 2 = X 2 -XI = \Y1 H ^ ~ a °)^' - a ?) ( 1T ) 

where Xo is the value of the minimal \ 2 characterized by the free fit parameters a°. In 
our fits we have d = 21 and x 2 is calculated by adding the total systematic and statistical 
experimental errors in quadrature. The uncertainties Aa^ = a,i — a® are constrained by 

A X 2 < T 2 (18) 

where the tolerance parameter T was chosen to be [57] 

T 2 = T 2 a = v / 2iV/(1.65) 2 , (19) 
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i.e. T ~ 4.7 since N = 1739 is the total number of data points considered in our global 
fits. The inversion involved in evaluating Ac^ in (17), subject to the constraint (18), 
was performed with the help of the normalized eigenvectors [58] of whose iterative 
calculation followed [59]. The calculation of all the uncertainties presented in our paper 
was performed according to the master equation (24) of [SH] whose particular implication 
for Acii is specified in eq. (30) of [58]. Our choice for the displacement distance t entering 
these latter equations was t = T, an assumption made in most subsequent publications 
and analyses. (When comparing our uncertainty results with the ones of CTEQ [21 [58] 
where T = 10 has been assumed, we rescale these CTEQ uncertainties according to 
Aci; — > 0.47Acij in order to comply with our T = 4.7 in (19).) 

As suggested in [2], we included in our final error analysis only those parameters 
that are actually sensitive to the input data set chosen, i.e. those parameters which are 
not close to 'flat' directions in the overall parameter space. With current data, and our 
functional form (5), 13 such parameters, including a s , are identified and are included in 
our final error analysis; the remaining ill-determined eight polynomical parameters Af 
and Bf, with uncertainties of more than 50%, were held fixed. 

3 Quantitative results and very small- £ predictions 

A representative comparison of our dynamical LO and NLO(MS) results with the relevant 
HERA(H1, ZEUS) data on the proton structure function F^x^Q 2 ) is presented in Figs. 
1 and 2. Due to our valence-like input, the small-x results (x £ 10~ 2 ) are predictions 
being entirely generated by the QCD Q 2 -evolutions. This is in contrast to a 'standard 
fit' where the gluon and sea input distributions in (5) do not vanish as x — > {a g ^ £j 0) 
at Qq = 2 GeV 2 . For comparison we have also performed such a standard fit shown by 
the dashed-dotted curves in Figs. 1 and 2. In both cases the data in Figs. 1 and 2 are 
well described throughout the whole medium- to small-x region for Q 2 >, 2 GeV 2 and 
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thus perturbative QCD is here fully operative. At Q 2 < 2 GeV 2 the theoretical results 
fall below the data in the very small-x region; this is not unexpected for perturbative 
leading twist-2 results, since nonperturbative (higher twist) contributions to F 2 (x,Q 2 ) 
will eventually become relevant, even dominant, for decreasing values of Q 2 . It should 
be emphasized that all of our valence-like input distributions at /i| Q = 0.3 GeV 2 and 
/^nlo = 0-5 GeV 2 as well as the ones for the 'standard fit' at Ql = 2 GeV 2 are manifestly 
positive. This is in contrast to negative gluon distributions in the small-x region observed 
in other standard fits [21 IH [5] • Furthermore the more restrictive ansatz of the valence-like 
input distributions at small-x as well as the sizeably larger evolution distance (starting at 
Qo < 1 GeV) imply smaller uncertainties concerning the behavior of structure functions 
in the small-x region than the corresponding results obtained from the common 'standard 
fits', in particular as Q 2 increases. This is illustrated in Fig. 3 for the NLO(MS) results 
where the error bands correspond to a la uncertainty. Since our valence-like sea input has 
a rather small power of x, i.e. vanishes only slowly as x — > 0, the uncertainties of the sea 
dominated F 2 (x, Q 2 ) turn out to be not too different from the standard fit where the sea 
increases as x — > (negative power of x) already at the input scale Qq = 2 GeV 2 . Notice 
that the uncertainties generally decrease as Q 2 increases due to the QCD Q 2 -evolutions 

[53 EH]. 

Our NLO(MS) valence-like input distributions at Ql = /i^Lo = 

0.5 GeV 2 together 

with their la uncertainties are shown in Fig. 4. They have been parametrized according 

to (5) with the parameters given in Table lj^] For comparison the GRV98 input [8] is 

displayed as well, which turns out to be very similar except for the gluon which peaks 

at a slightly larger value of x. However, such differences are merely within a 2a band 

of our new results. The valence-like gluon input at low Q 2 < 1 GeV 2 in Fig. 4 implies 

2 It should be mentioned that there is a correlation between the (chosen) value of a s (M|) and the 
resulting values of the valence-like input scales /ilo,nlo which increase with a s (Af|). Since we did not 
want to fix a s (M|) at the LEP value of 0.118, we performed fits for various fixed values of n by imposing 
a valence-like input structure (a g u+t i > 0) and keeping a s as a free fit parameter. Then we fixed the 
best choice for /z (/i 2 = 0-3 GeV 2 , I^lo = 0-5 GeV 2 ) and performed the final precision fits and error 
analyses. 



9 



a far stronger constrained gluon distribution at larger values of Q 2 as compared to a 
gluon density obtained from a 'standard fit' with a conventional non- valence-like input 
at Q 2 > 1 GeV 2 as can be seen in Fig. 5. As already mentioned, this is in contrast to 
the sea distribution u + d in Fig. 5 where the valence-like sea input in Fig. 4 vanishes 
very slowly as x — > and thus is similarly increasing with decreasing x down to x ~ 0.01 
as the sea input obtained by a standard fit. Therefore the la uncertainty band of our 
dynamically predicted sea distributions at larger values of Q 2 in Fig. 5 is only marginally 
smaller than the corresponding one of the standard fit. The relevant input parameters of 
our 'standard fit' can be found in Table 2. As expected for the dynamical fit, starting 
from a low input scale with valence-like distributions, a s (M|) in Table 1 is somewhat 
stronger constrained due to the larger evolution distance than the corresponding result of 
the standard NLO(MS) fit in Table 2. Keeping in mind that our stated errors always refer 
to lcr uncertainties, our standard fit error of 0.0021 for a s (M|) in Table 2 is compatible 
with the 2a uncertainty also stated in the literature (see, e.g., [2] and the discussion in 
[3]). It should be furthermore mentioned that our NLO(MS) results for a s (M|) in Table 
1 and 2 are, within about a la uncertainty, also compatible with the ones obtained from 
analyzing only DIS structure functions (for a recent summary, see |60j). 

At this point it should be mentioned that the standard CTEQ [2] fit resulted, very 
surprisingly, in a valence-like input gluon distribution at a scale as large as Qq = m 2 ~ 1.7 
GeV 2 . Thus this CTEQ6 gluon distribution is expected to be similarly tightly constrained 
at Q 2 > Qq as our dynamical results starting from a valence-like input at Qq = 0.5 GeV 2 . 
That this is indeed the case is illustrated in Fig. 6 where the la uncertainties of the CTEQ6 
gluon [2J are similar in size to our dynamical results, whereas a common 'standard fit' 
(being based on an increasing input distribution as x — > 0) results in a sizeably larger 
uncertainty. The situation is different for the sea distribution in the small-x region; here 
the CTEQ6 input at Qq = m 2 increases at x — > as expected, and thus the la uncertainty 
is comparable to our 'standard fit' result as shown in Fig. 7 - both being larger than the 
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la uncertainty obtained from our dynamical fit based on valence-like inputs at Qq = 0.5 
GeV 2 . 

The heavy charm and bottom quark contributions to F 2 at LO and NLO(MS) are 
compared with recent HERA data in Figs. 8 and 9. The impressive agreement with 
present measurements for F 2 C and _F 2 b illustrates that the n/ = 3 FFNS is entirely reliable. 
As already discussed in Sec. 2, the NLO results are rather insensitive to the chosen 
factorization scale fip (^p = 4m^ or fi 2 F = Q 2 + 4m^). Again the la uncertainties of 
these dynamical predictions are distinctly smaller than the ones implied by the standard 
fit. It should be furthermore reemphasized that within the FFNS heavy quarks (h = 
c, b, t) are always produced as final states in fixed order perturbation theory via hard 
production processes initiated by the light partons of the nucleon (u, d, s quarks and 
the gluon g). The perturbative stability of heavy quark production (S3] as well as the 
agreement with experiment in Figs. 8 and 9 even at Q 2 3> m\ indicate that there is no 
need to resum supposedly 'large logarithms' (lnQ 2 /m^), which is of course in contrast to 
genuine collinear logarithms appearing in light (massless) quark and gluon hard scattering 
processes. Therefore only the nf = 3 light u,d,s quark flavors and gluons constitute the 
'intrinsic' genuine partons of the proton and the heavy c, b, t quark flavors should not be 
included in the parton structure of the nucleon, not even at Q 2 3> ml [53]. However, 
somewhat dissenting views were recently summarized in [61] . 

The measurements of Drell-Yan dilepton production in pp and pd collisions [HJ H2] 
are instrumental in fixing A = d — u (or d/u) [62]. In Fig. 10 we display our dynamical 
NLO(MS) result for a pd /2a pp together with the ±la uncertainty band as well as the previ- 
ous GRV98 result which agree in the statistically relevant x-region, with x 2 referring to the 
average fractional momentum of the target partons. Note that a pN = d 2 a pN /dxpdM 2 +fi _ 
with xp = X\ — x 2 . In LO a pN ~ ^ urfs e^ [9(^1)^(^2) + 9(^2)9(^1)] where X\ and x 2 refer 
to the Bjorken-x of the quarks in the beam (p) and nucleon target (N), respectively. Ex- 
perimentally xp > (xi > x 2 ) and consequently the Drell-Yan cross section is dominated 
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by the annihilation of a beam quark with a target antiquark. For x\ ^> x 2 one obtains 
a pd/ 2( jPP ~ 1 [1 + d(x 2 )/u(x 2 )] at a scale Q 2 = M 2 + ^_ in g(x 2 , Q 2 ). 

Finally the pp Tevatron high-px (or E^) inclusive jet data HZ] are compared in 
Fig. 11 with our dynamical LO and NLO(MS) results, as well as with the ones of CTEQ6 
[2]. The small lex error bands are almost invisible on the huge logarithmic scale used. Our 
NLO result almost coincides with the one of CTEQ. There is a clear improvement at NLO 
as compared to LO which falls slightly below the data at pr ~ 300 GeV. Nevertheless the 
LO high-px fit corresponds to x 2 /dof ~ 1 which is only twice as large as at NLO. 

As discussed in Sec. 2 we have explicitly used for our analysis the experimentally 
directly measured 'reduced' DIS cross sections (7) which, for not too large values of Q 2 , 
are dominated by the one-photon exchange cross section a r = F 2 — (y 2 /Y + )Fi J where 
y = Q 2 /xs. The importance of using this quantity has recently been emphasized [63] 
since the effect of becomes increasingly relevant as x decreases at a given Q 2 where y 
increases. This is seen in the data as a flattening of the growth of a r (x, Q 2 ) as x decreases 
to very small values, at fixed Q 2 , leading eventually to a turnover (cf. Fig. 12). At the 
lower values of Q 2 in Fig. 12 it was not possible in [63J to reproduce this turnover at 
NLO. This was mainly due to the negative longitudinal cross section (negative F^(x, Q 2 )) 
encountered in [63] . Since all of our cross sections and structure functions are manifestly 
positive throughout the whole kinematic region considered, our dynamical NLO (MS) 
results in Fig. 12 are in good agreement with all small-x HERA measurements [2B1I23J- For 
completeness we compare in Fig. 13 our dynamical (leading twist) NLO (MS) predictions 
for Fl(x,Q 2 ) with a representative selection of (partly prelinimary) HI data [2E1 EI] at 
fixed W ~ 276 GeV. The standard fit result with its sizeably larger ±1<t error band is, 
for comparison, shown as well. Our NLO results for Fl, being gluon dominated in the 
small-x region, are in full agreement with present measurements, which is in contrast 
to expectations [31 |63] based on negative parton distributions and structure functions 
at small x. To illustrate the manifest positive definiteness of our dynamically generated 
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structure functions we show F^(x, Q 2 ) in Fig. 13 down to Q 2 = 1 GeV 2 although a leading 
twist-2 prediction should not be confronted with data below, say, 2 GeV 2 . 

As our parameter-free small-x predictions for parton distributions at x < 1CT 2 are 
entirely of Q CD-dynamical origin and depend rather little on the detailed input param- 
eters at x 10~ 2 , it is interesting to study these predictions in kinematic regions not 
accessible by present DIS experiments. Of particular interest are, as already emphasized 
in the Introduction, calculations of weak ^uN cross sections of ultrahigh energy cosmic 
neutrinos [12, 15, 65-69] which afford a (reliable) knowledge of parton distributions at 
the weak scale Q 2 = down to x ~ 1CT 9 for highest energies E v ~ 10 12 GeV. This 
requires extrapolations into the unmeasured small- x region x < 10 -3 . Since F%{x, Q 2 ) is, 
in the very small-x region, dominated by u(x, Q 2 ) and d(x, Q 2 ), as are the CC neutrino- 
(isoscalar) nucleon cross sections, the -Ff structure function can be utilized for estimating 
the magnitude of uncertainties of the predictions in the extreme small- x region which 
are shown in Fig. 14. At Q 2 = our dynamical NLO predictions correspond to a 
±1(T uncertainty of about ±7% at x — 10~ 9 whereas the uncertainty of the extrapolation 
of a standard fit is about twice as large. At smaller scales the uncertainties obviously 
increase as illustrated in Fig. 14 at Q 2 = 100 GeV 2 . Taking into account previous ex- 
trapolation ambiguities j8], one can conclude [12] that the dynamically predicted small-x 
parton distributions allow neutrino-nucleon cross sections to be calculable with an accu- 
racy of about 10% at highest cosmic neutrino energies. It should be mentioned that an 
ad hoc fixed power law of x extrapolation of the standard CTEQ6.5 structure functions 
[70] to x — 10~ 8 at Q 2 = [7T] lies, accidentally, only about 10% below our dynam- 
ical NLO prediction in Fig. 14. On the other hand, an alternative parametrization [71] 
of present HERA(ZEUS) data which is not QCD oriented but based on analyticity and 
unitarity gives, when extrapolated to x — 10~ 8 , a factor of about 6 smaller a value for 
F% (10~~ 8 , Myy) than our prediction in Fig. 14. Since the perturbative dynamical QCD 
predictions for the small-x behavior of structure functions down to x — 10~ 5 proved to 



13 



be in agreement with later HERA measurements as discussed in the Introduction, it is 
hard to imagine that perturbative QCD dynamics and evolutions should become entirely 
inappropriate at x — to 10" 9 at even much larger scales. 

In order to test the dependence of our results on the specific choice of the factorization 
scheme in NLO, we have repeated our dynamical analysis in the DIS factorization scheme 
as outlined in Sec. 2. Since the (light) parton distributions in the DIS scheme are defined 
via the F 2 structure function in (12), it is not very surprising that, in contrast to the DIS 
parton distributions themselves [72], the results for physical observables directly related 
to DIS structure functions are very similar. Indeed the results in the DIS factorization 
scheme for F^x^Q 2 ) in Figs. 1 to 3 and 14 are practically indistinguishable from the 
ones in the MS scheme, as are the results for F^^x, Q 2 ) in Figs. 8 and 9 and the ones 
for <j r (x,Q 2 ) and F^x^Q 2 ) in Figs. 12 and 13, respectively. Differences become visible 
only for processes which are not directly related to DIS structure functions such as Drell- 
Yan cross sections but there the differences lie well within the ±lcx uncertainty of the 
NLO(MS) results as illustrated in Fig. 10 for the DY-asymmetry. The parameters for the 
input parton distributions in the DIS scheme and the corresponding a s (M|) can be found 
in Table 1. It should be furthermore emphasized that our results and predictions are also 
stable to within less than about 20% when compared to previous analyses and fits. This is 
illustrated in Fig. 15 for our present dynamical NLO (MS) results when compared with our 
previous GRV98 results [8] . The situation is similar for more recent and previous standard 
CTEQ and MRST parton distributions for their relevant ranges in x, and holds also by 
comparing CTEQ and MRST distributions with each other [21 El E7J EU E3J EU ■ It should 
be emphasized that heavy quark mass effects have always been fully taken into account 
in our previous [7J [S] and present analyses. This is in contrast to the previous CTEQ6 
analysis [2] where charm has been treated in the zero-mass approximation. The recent 
inclusion of finite charm mass effects in CTEQ6.5 [70] reduces the charm contribution to 
F 2 (x,Q 2 ) which is compensated by larger u = u v + u and d = d v + d distributions at 
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small x as compared to CTEQ6 [2J. That such an 'enhancement' has always been present 
in our dynamical u and d distributions is illustrated in Fig. 16, since our present and 
previous (cf. Fig. 15) distributions differ very little from the CTEQ6.5 ones. Therefore 
our predicted hadronic W ± /Z° production cross sections, for example, at Tevatron and 
LHC are similar to the 'enhanced' ones observed in [70] . 

Of course more recent parton distributions have a higher precision due to the higher 
statistics of the data, but we have not experienced essential qualitative and quantitative 
changes during the past decade. It is reassuring to see that our knowledge of the fun- 
damental partonic structure of matter has essentially remained unchanged over the past 
years. 



4 Summary and Conclusions 



Utilizing recent DIS measurements and data on Drell-Yan dilepton and high-.©! inclusive 
jet production, we have redone a previous [5] global fit for the dynamical parton distribu- 
tions of the nucleon in the LO and NLO of perturbative QCD. The small-x (x £ 10 -2 ) 
structure of dynamical parton distributions is generated entirely radiatively from valence- 
like, manifestly positive, input distributions at an optimally chosen input scale Qq < 1 
GeV. The NLO results are stable with respect to a different choice of the factorization 
scheme (MS versus DIS). The predictions for the longitudinal structure function F L (x, Q 2 ) 
at small x, for example, are positive throughout the whole kinematic region considered, 
in agreement with (partly preliminary) data. We have augmented our analyses with an 
appropriate uncertainty analysis and found that the newly determined dynamical distri- 
butions are compatible with the former [7J IE] ones, where heavy quark mass effects have 
always been fully taken into account. The stability of these results guarantees a reliable 
calculation of cross sections for, e.g., heavy quark, W ± , Z°, and high-pr jet production 
at hadron colliders like Tevatron and in particular LHC. 
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Our dynamical distributions have also been compared with conventional ('standard') 
ones obtained from non-valence-like positive input distributions at some arbitrarily cho- 
sen higher input scale Qo > 1 GeV. For this purpose we have performed a 'standard 
fit' as well, assuming Qq = 2 GeV 2 . The uncertainties of these latter distributions are, 
as expected, larger, in particular in the present (experimentally) unexplored extremely 
small-x region relevant for evaluating ultrahigh energy neutrino-nucleon cross sections in 
astrophysical applications. Here we provide predictions down to x ~ 1CT 9 at the weak 
scale Q 2 = as required [121 E3 [66] for highest cosmic neutrino energies of 10 12 GeV. 
These predictions are strongly constrained within the dynamical parton model and are 
entirely of QCD-dynamical origin in the very small-x region. Furthermore, as mentioned 
in the Introduction, previous predictions [El |9] for the small-x region based on the dy- 
namical parton model and the data available at that time were subsequently confirmed 
[TOl fTT] at HERA. The presently available very precise small-x data [28], [29J utilized here 
allows us to be quite confident about the reliability of our improved small-x predictions 
within the framework of the successful dynamical parton model. 

A FORTRAN package (grid) containing our new dynamical LO, NLO(MS), and 
NLO(DIS) parton densities, the light (u,d,s;g) F2 Sh± (x,Q 2 ) as well as F 2 c,b (x, Q 2 ), cal- 
culated in the fixed order FFNS, can be obtained by electronic mail or on request. The 
NLO(MS) uncertainty estimates will be also included. 
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Table 2: As Table 1 but for 
input scale Q\ = 1 GeV 2 . 
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Figure 1: Comparison of our dynamical (dyn) LO and NLO(MS) as well as standard 
(std) NLO small-ac results for F$(x,Q 2 ) with HERA data for Q 2 > 1.5 GeV 2 [2E1 122] ■ 
The parameters of the valence-like input distributions for the dynamical predictions are 
given in Table 1 and the ones for the standard results in Table 2. To ease the graphical 
presentation we have plotted F$(x, Q 2 ) + 0.5 x i{Q 2 ) with i{Q 2 ) indicated in parentheses 
in the figure for each fixed value of Q 2 . 
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Figure 3: Typical ±1<t uncertainty bands of our dynamical and standard NLO(MS) results 
in Fig. 1 for two representative values of Q 2 . To ease the visibility of the two error bands 
at Q 2 = 2.5 GeV 2 we have subtracted 0.2 from the stdNLO result as indicated. For 
illustration two HI and ZEUS data points from Fig. 1 with their almost invisible errors 
are shown as well at Q 2 = 2.5 and 15 GeV 2 , respectively. 




Figure 4: The valence-like input densities together with their ±ler uncertainties at = 
/^nlo = 0-5 GeV 2 for our dynamical NLO(MS) results. The central curves follow from (5) 
with the parameters given in Table 1. The strange sea s = s vanishes at the input scale. 
The NLO GRV98 input [8] is also shown by the dashed curves for comparison. 




Figure 5: The small-x NLO(MS) predictions of our dynamically (radiatively) generated 
gluon and sea-quark distributions, together with their ±lcr uncertainties, as compared 
to the results of a standard fit. To ease the visibility of the two error bands of x(u + d) 
we have multiplied the stdNLO results by 0.9 as indicated. The corresponding GRV98 
predictions [8] lie within the la band of our new dynNLO results. 
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Figure 6: Comparing the ilex error bands of our dynamical (dyn), standard (std) and 
CTEQ [2] NLO(MS) gluon distribution at small x for various fixed values of Q 2 . Note 
that Q 2 = 2 GeV 2 is the input scale of the standard fit which is close to the CTEQ input 
scale Qq = m 2 ~ 1.7 GeV 2 , where the standard CTEQ6 fit employs a valence-like gluon 



input (i.e., xg(x,Ql) — > as x — > 0). Due to the sizeably different input scales, the 
CTEQ6 gluon falls up to 30-40% below our dynNLO gluon for x < 10" 3 and Q 2 > 10 
GeV 2 . The results at Q 2 = 2 and 20 GeV 2 have been multiplied by 0.5 and 2, respectively, 
as indicated in the figure. 
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Figure 7: As in Fig. 6 but for the sea-quark distribution x(u + d). Notice that here the 
CTEQ input distribution at Ql ~ 1.7 GeV 2 is not valence-like (i.e., x{u + d)(x,Ql) 
for x — > 0). 




Figure 8: The dynamical NLO(MS) predictions for F% in the strict n/ = 3 FFNS, choosing 
fi 2 F = Am 2 , with m c = 1.3 GeV, together with the ±1<t uncertainty band. For comparison 
we also display the central LO predictions which are entirely due to the 7*-gluon fusion 
subprocess 7*g — > cc. The charm production data as obtained from D* measurements are 
taken from EH] (solid and open squares) and the HI direct track measurements from 
PH] (open circles). 
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Figure 9: As in Fig. 8 but for F$ with m& = 4.2 GeV and the bottom production data 
taken from [40J. 
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Figure 10: Our dynamical NLO result in the MS factorization scheme, together with 
its ±lcr uncertainty, for a pd /2a pp appearing in the Drell-Yan asymmetry Aoy = {v pp — 
a pn ) I (a pp +a pn ) as a function of the average fractional momentum x<i of the target partons. 
The GRV98 NLO MS-result [5] is shown for comparison. The dynamical NLO(DIS) result 
in the DIS factorization scheme is shown by the dashed-dotted curve. The data for the 
dimuon mass range 4.6 < M^+^- < 12.9 GeV are taken from [32] ■ 




Figure 11: The pp Tevatron high-pp inclusive jet data HTj compared with our 
namical LO and NLO(MS)results, as well as with the NLO CTEQ6 result [2]. 




Figure 12: The dynamical NLO(MS) predictions, together with their ±lcr uncertainties, 
for the 'reduced' DIS cross section a r (x, Q 2 ) = F 2 — (y 2 /Y + )F L . The HERA data for some 
representative fixed values of Q 2 are from [281 EH] . 




Figure 13: Dynamical and standard NLO(MS) results for F^(x,Q 2 ) together with their 
±ler uncertainty bands. The (partly preliminary) HI data [281 El] are at fixed W ~ 276 
GeV. 
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Figure 14: Dynamical NLO(MS) predictions for F%(x, Q 2 ), together with their la uncer- 
tainties, for extremely small values of x. The la uncertainties of the standard (std) NLO 
fit extrapolations are shown by the vertical bars. The dotted curves are the contributions 
from the light (n/ = 3) quarks and gluons to Ff for the dynamical (dyn) NLO result. In 
other words, the difference between the dotted and solid curves is due to NLO heavy quark 
(charm, bottom) contributions which derive from photon-gluon (quark) fusion processes. 
The dynamical GRV98 predictions [8] lie within the ±la band of our present dynNLO 
predictions. 




Figure 15: Comparing the present dynamical dynNLO(MS) parton distributions with the 
previous ones of GRV98 [g] at Q 2 = 10 GeV 2 . 
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Figure 16: Comparing the present dynamical NLO(MS) u = u v + u and d = d v + d 
distributions with the ones of CTEQ6.5 [70J at Q 2 = 10 GeV 2 . These ratios remain 
practically unchanged at higher scales, like Q 2 = relevant for W ± production. The 
shaded areas represent the estimated ±lcr uncertainty band of our dynNLO analysis. 
Notice that in the relevant small-x region these ratios would be practically unaltered if 
the GRV98 distributions [8] were used instead of the dynNLO ones, since dynNLO/GRV98 
~ 1 as evident from Fig. 15. 



